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MINIMAL POLYNOMIALS FOR THE COORDINATES OF THE 

HARBORTH GRAPH 

EBERHARD H.-A. GERBRACHT 



Abstract. The Harborth graph is the smallest known example of a 4-regular 
planar unit-distance graph. In this paper we give an analytical description 
of the coordinates of its vertices for a particular embedding in the Euclidean 
plane. More precisely, we show how to calculate the minimal polynomials of 
the coordinates of its vertices with the help of a computer algebra system, and 
list these. Furthermore some algebraic properties of these polynomials, and 
consequences to the structure of the Harborth graph are determined. 



o 
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, 1. Introduction 

The Harborth graph (see Figure 1) is the smaUest known example of a 4-regular 
planar unit-distance graph. That is a planar graph, all of which edges are of unit 
^■f-^ . length, with exactly four edges meeting in each vertex. This graph was named after 

^ ' its discoverer H. Harborth, who first presented it to the general public as a research 

i problem in ^4^ and to a large international audience in a talk at the Eugene Strens 

' Memorial Conference on Recreational Mathematics and its History in 1986 [5l I12j. 

At both occasions he posed the question, if a smaller example of a 4-regular planar 
' unit-distance graph could be found. 

\^ . Curiously enough, up until now nearly all published pictures of the Harborth 

' graph - even the original ones in [H [5] , as well as those in textbooks [H [7] - seem 

fH , to be slightly vague and inaccurate, with the vertices always being depicted by 

large dots. Furthermore there has not been given any analytical description of 

a"* . the Harborth graph, yet. That is, if we consider the graph as embedded in the 

Euclidean plane with a given coordinate system, the coordinates of the vertices 
^ . have never been calculated exactly. This has gone to the point that the Harborth 

graph was even thought by some to be nonrigicU, which, as the results of this paper 
imply, cannot be the case. Thus the wish for an "exact" description remained. This 
' wish has recently been expressed again in the world wide web , which prompted 

this paper. 

With the advent of dynamical geometry systems, several authors were finally 
able to produce more precise pictures [Ml [11] , which led to further evidence that 
there is one unique realisation of the Harborth graph in the Euclidean plane. With 
this paper we go one step further: using one particular way of construction, we 
set up a set of quadratic equations which completely describe the coordinates of 
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^ If we consider the Harborth graph as a (mechanical) framework consisting of rigid bars 
interconnected by rotable joints, nonrigidity means that some vertices can be moved with respect 
to each other, that is the whole framework allows motions which are different from congruences 
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Figure 1. The Harborth graph embedded in the Euchdean plane. 



crucial vertices of the Harborth graph. Using this initial set of equations, we will 

show that all the coordinates arc algebraic numbers, and wc will calculate their 
minimal polynomials with the help of a computer algebra system. Even though, as 
we will see later on, it is impossible to solve the corresponding algebraic equations 
exactly, which in our understanding means in terms of radicals, nevertheless wc thus 
have achieved an exact analytic description of the Harborth graph, since, together 
with easy to calculate numerical approximations for the actual coordinates, these 
polynomials uniquely determine each coordinate. 

The author is indebted to C. Adelmann and H. Lowe from the Technical Uni- 
versity Braunschweig, Germany, for a number of very valuable discussions on the 
subject. Furthermore he owes thanks to the Institute for Mathematical Physics of 
the TU Braunschweig, especially to the research group of R.F. Werner, for making 
its computing facilities available to him for this research. 

2. Using Dynamic Geometry Software and Numerical Analysis 

2.1. Geometric construction. Because of the obvious twofold symmetry of the 
Harborth graph, it is enough to analyse one of its quarters. Therefore, in a first 
step we construct one of these quarters (see Figure 2), using one of the existing 
(imprecise) first generation pictures as a blueprint. 

We start from an initial isosceles triangle ABC of fixed but arbitrary height 
T, with two sides being of unit length, and a neighbouring symmetric trapezoid 
BCDE, which has the side EC with the initial triangle in common. The parallel 
sides of the trapezoid are chosen to be of length 2 and 3 respectively. The remaining 
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points are constructed from this initial configuration by using compass and ruler 
techniques. In the following, we list the necessary steps. Thereby Circ(P, r) denotes 
the circle with center P and radius r, and H the operation of letting two geometric 
figures intersect. Thus we get 



Although the intersection of two circles usually consists of two points, we use this 
notation as if there was no ambiguity. We are allowed to this, because we choose the 
resulting points of intersection according to our blueprint. Thus, e.g., F is chosen 
in such a way that the quadrangle ABEF is convex. In the sequel we will call the 
configuration of the points A to J thus constructed the Harborth configuration. 
We will do this, even if the parameter T initially has not been chosen correctly, so 
that the configuration cannot be completed to the whole Harborth graph. 

Clearly, proceeding as described, with using an arbitrary nonnegative value for 
the height T of the initial triangle, the line through the final crucial points H and 
J will only by (a very small) chance meet the line through A and C at an angle 
(f of 90°. But, in order that the Harborth configuration can be completed to the 
whole Harborth graph, we have to make sure that this last condition is satisfied. 

Using the dynamic geometry software GeoGebra [8], we are able to manipulate 
the Harborth configuration in dependence on the parameter T. Furthermore with 
GeoGebra we are able to read off approximate values (up to five decimal places) 
for the resulting angles and coordinates. At this early stage of our investigations 
already, some important observations can be made: 

Remark 2.1. 

(1) The point F, and thus the whole Harborth configuration only exists for 

T e [0, 6], where h is approximately 0.13504. 

(2) For T e [0,6], the angle tp between AC and HJ lies approximately in the 

interval [85.88496°, 94.59043°]. 

(3) For T ~ 0.12073, the angle <p approximates 90°. 

The upper bound b in the above remark is determined by the fact that for T > b 
the distance between A and E becomes greater than 2, and the circles with centers 
A and E of radius 1 do not meet anymore. 

As hinted at in the introduction, once we have calculated the exact values for the 
extremal angles - which we will do in one of the following sections - the observations 
collected in Remark |2 . 1 1 allow us to give the first analytic prooQ of the planarity of 
the Harborth graph which does not resort to pictures or models only. 

Theorem 2.2. There exists T G [0,6] such that (p — 90° exactly, and the Harborth 
configuration can be completed to the Harborth graph. 

Proof. lOj Since the geometric operations we used in our construction (i.e., drawing 
circles and letting them intersect) depend continuously on their parameters (i.e.. 



The basic idea of this proof was first communicated to the author by H. Lowe in 2003. Here 
it is presented with his kind permission. 



(2.1) 
(2.2) 
(2.3) 
(2.4) 
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Circ(A, 1) n Circ(£;, 1) 
Circ(i^, 1) n Circ(i:>, 2) 
Circ(L>,2) n Circ(G,2) 
Circ(i^, 1) n Circ(G, 1) 
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Figure 2. A quarter of the Harborth graph, created with GeoGe- 
bra [5] - the point B (and thus the height T) chosen in order that 
ip = 90°, i.e., the configuration can be completed to the whole 
Harborth graph. 



centers and radii), and since the composition of continuous functions is again con- 
tinuous, the angle ip depends continuously on the height T as long as the Harborth 
configuration exists. By the Intermediate Value Theorem the existence of at least 
one T in the above interval [0, b] is assured such that p is precisely 90°. □ 

Corollary 2.3. The Harborth graph indeed is a 4-regular unit-distance graph which 
is planar. 

From now on, our main goal will be to determine a precise description of this 
particular T, the existence of which we have shown above, without resorting to 
trial-and-error. 

2.2. Setting up algebraic equations. To describe the points of the Harborth 
graph more precisely we need to introduce coordinates. For the time being, we 
choose the point A as the center of our coordinate system, and the ray which extends 
the base side AC of the triangle ABC as the positive x-axis. When considering the 
Harborth graph as a whole, this is by far not the most obvious choice. In fact, in 
a later section we will use its center of symmetry as the origin of a more natural 
coordinate system, and then will have to "translate" our intermediate results. As 
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we proceed, we will see that in some places this will prove to be quite cumbersome, 
which will retroactively justify our initial choice of coordinates. 

Now let {t,T) be the coordinates of the point B. Since A, B, and C form an 
isosceles triangle, the coordinates of the point C are given by (2t, 0). 

Denoting the coordinates of any point P of the graph different from B with 
xp, Up), we have: 



(2.5) 
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{xe - xpf + (yE - ypf 
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{xF - xaf + {yp - ycf 


- 1 
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(2.13) 


{xD - Xaf + {yo - ycf 


-4 
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(2.14) 


{xD - xnf + [yo - ynf 


-4 
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(xG - XHf + {yo - yuf 


- 4 
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(2.16) 


{xF - xjf + {yp - yjf 


- 1 
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(2.17) 


{xG - xjf + [yc - yjf 


- 1 


= 


(2.18) 


XH - 


xj 


= 



Let us shortly comment on the meaning of these equations: (|2.6p . (|2.7p . respec- 
tively (|2.8p , (|2.9p , define the vertices D and E of the trapezoid which do not belong 
to the initial triangle ABC. Equations (|2.6p and (|2.8p stem from the fact that the 
line BC meets CD at an angle of 60° and the hnc EB meets BC at 120°. Every 
further pair of Equations (|2.10p - (|2.17p is chosen in accordance to the geometric 
constructions described in (|2.ip - (|2.4p . each pair defining one of the points F - J. 
Finally, Equation (|2.18p has to be satisfied in order that the lines given by AC and 
HJ meet at an angle of 90°. 

Using this set of equations, and approximations for the coordinates, which we 
read off from Figure 2 with the help of GeoGebra, we are already able to calculate 
arbitrary precise approximations of all coordinates by using standard numerical 
algorithms. E.g., the results below show the coordinates with an exactness of 15 
digits; they were calculated with Mathematica, version 4.0.1.0. 

B = {t,T) ~ (0.992685948824186,0.120725337054926) 

D ~ (0.809996600722107, 2.760161754567202) 

E ~ (0.209102417540010, 1.960833173433061) 

(2.19) F ~ (-0.061398137844065,0.998113354619244) 

C ~ (-0.838419516770942, 1.627587561152422) 

H ~ (-0.995049481192288, 3.621444891616507) 

J ~ (-0.995049481192288, 0.639930204451542) 
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3. Deducing an equation for the ?/-coordinate of the point B 

In this section we will deduce the minimal polynomial for the coordinate ys ~ 
T, i.e., the unique primitive integer polynomial Pt of smallest degree such that 
Pt{T) = holds [2, Section 4.1.1]. 

Let us shortly describe the main approach which we will repeatedly take: Given 
two polynomial equations, which are satisfied by certain coordinates, we will cal- 
culate the resultant of the corresponding polynomials with respect to one of the 
appearing variables, thus eliminating this particular variable. Sometimes, if the 
polynomials are not too complicated, we will use Groebner basis techniques to 
treat more than two polynomial equations simultaneously. Both procedures lead 
to new polynomials or polynomial equations, which are consequences of the orig- 
inal one^, but contain fewer variables. To keep expressions from becoming too 
complicated, and running times from becoming too long we will try to factor the 
resulting polynomials. Many times it will prove to be advantageous to allow the 
factorization to be done over the ring extension Z[V3]. If a particular polynomial 
is reducible we will continue our deliberations with that factor, which corresponds 
to the actual values of the coordinates. To check this, we use numerical approxi- 
mations analogous to those given in Section [2.21 but which are precise to an error 
of e = 10"^"°. Most calculations, especially those of resultants, factorizations, and 
numerical evaluations, were done with Mathematica, version 4.0.1.0. 

The succession of eliminations will be determined by the order in which the cor- 
responding points were constructed. E.g., in the initial isosceles triangle ABC due 
to the choice of the coordinate system all point coordinates are directly expressible 
in terms of T. Next we will determine polynomials which describe the connection 
between the coordinates of D and E, respectively, and the parameter T . After 
that the polynomials for the coordinates of F are calculated by using those for the 
coordinates of D and E and eliminating the variables in between. In principle, 
continuing this procedure would lead to polynomials in xh and T, respectively xj 
and T. Using the final equation xh — x,j, one should be able to deduce one poly- 
nomial in the variable T alone. Unfortunately, due to the increasing complexity 
of expressions we were not able to continue this line of thought to its conclusion, 
but had to resort to an alternative way. Nevertheless wc will try to push as far as 
possible with this approach, and come up with an alternative, when it proves to be 
necessary. 

In the sequel, we will switch between listing the polynomials and the correspond- 
ing polynomial equations at will. When only a polynomial P is given it should be 
understood that the coordinates appearing in P satisfy the corresponding polyno- 
mial equation P = 0. 

3.1. Prom A to F. Using the procedure GroebnerBasis with Equations (|2.5|) - 
((27)) as input, like 

GroebnerBasis [{t ~2+T~2-l , -t* (xD-2*t) +T*yD-3/2 , (xD-2*t) ~2+yD~2-9} , 
{t,xD,yD,T}] 

one of the polynomials we get is 

(3.1) Py^^T 27 - 36T2 + 12T • ijd - Ayl. 



^For a more precise description see, e.g., [3] or [9]. 
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Analogously, by changing the order of variables, 

0=1- 56T^ + 784r^ - 8x1 - 208T^ ■ x^ + 16x% 

can be deduced. This last equation is irreducible over Z, but factors over Z[V3] 
into polynomials, which are quadratic in xd '■ 

(3.2) (-1 + 28T2 - UVSTXD + 4x|,) • (-1 + 28T'^ + uVSTxd + 4a;|,). 

Using numerical results for T and xo in analogy to (|2.19|) . we see that only the 
first of these polynomials 

(3.3) P^^,T -1 + 28^2 - 12V3TXD +4x1, 

leads to the correct result. Solving p.3p and p.ip for xd and yo, respectively, 
and again discarding those solutions which do not describe the correct coordinates, 
we get explicit descriptions for the coordinates of D in terms of the parameter T: 

(3.4) XD = i (sVST + - ^2 

3 



(3.5) VD = ^[T + VsVl-T^ 

Starting with Equations (|2.5|) . ()2.8|) and ()2.9|1 and proceeding in the same manner 
as above, we are led to those polynomials which describe the dependence of the 
coordinates of the point E on the parameter T: 

(3.6) 3r2 - xl 
and 

(3.7) Py,,T:^-S + 7T^-4TyE + yl. 
Again, the corresponding equations can be explicitly solved: 

(3.8) XE = VST, 

(3.9) ye = 2T + Vs^l -^2. 

Next we continue by calculating the coordinates for the point F, once more using 
Mathematica's GroebnerBasis function. This time we start with the newly found 
set of Equations (|3?6| and ((37| . together with the Equations ((2T0| . (|2TT|1 defining 
F. From this we get the following polynomials which describe the dependence of 
the coordinates xf and yp of F on the parameter T : 

-81 + lOSOOT^ - 422496T* + 4272S84T^ - 19194112r** + 45801472T^'' - 

63111168T^^ + 48234496r^'' - 16777216r^^ + 
(1296 - 92448r^ + 1645056r'' - 9573888r^ + 30072832T* - 57655296ri° + 

(3.10) 66060288T^^ - 29360128r^'') x% + (-7776 + 228096T^ - 1555200T^ + 
5271552T^ - 12189696r** + 15728640T^'' - 9437184r^^) xp + 

(20736 - 152064T^ + 331776T'' - 196608T^ - 1310720r** + 2097152r^°) x% + 
(-20736 + 110592T^ - 442368T'' + 786432T'^ - 1048576T**) x%, 



(3.11) 



Py^ rj, 81 _ 648r^ + 144T'' - 2304r'^ + 4096r® + 
(432r - 864r^ + 6528T^ - 10240T'^) yp + 
(-216 + 1584T^ - 5376r'' + 9216r®) yj- + 
{-576T + 1536r^ - 4096r^) yf + (144 - 384T^ + 1024r'') yj.. 

Once again, the first of these polynomials factors over Z[-\/3] into two polynomials 
of total degree 8 and degree 4 in a;_F. Using the numerical values for T and xf, we 
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can deduce that only one of these describes the connection between the variables 
T and xf- It is 

Frp^T := - 9 + 600r^ - 3472r'' + SSSST*^ - 4096T* - 
(3.12) 8v^ (9r + 96r3 - 112T5 + 256r^) x^+S (9 - ISOT^ + 256r* - 640T'^) x%-+ 
8^ (36r - 96T^ + 256T^) x|. + 16 (-9 + 24r^ - x%. 

As said before, trying to continue like this to calculate polynomials for the remaining 
points G, H and J, will lead into a deadend, because the resulting equations become 
too unwieldy to handle, and take too much time to calculate, even with the help of 
Mathematica. Still, our main goal remains to find one single equation describing 
the parameter T alone. Consequently we have to take a step back, and use a slightly 
more indirect approach, which we will describe in the section following the next one. 

3.2. Interlude: Calculating the extremal values for which the Harborth 
configuration exists. With Equations p.Sp and p.9p thus available, we are able 
to calculate the exact maximal value for T, hinted at in Remark 1 2. II To this end, 
we first observe that in case of T being maximal the line segments AF and FE 
together, again form a straight line segment of twice the original length. Therefore, 
for T maximal, the coordinates of the point E satisfy x'^+y^ — A = 0. This, together 
with p.Sp and p.9p . after some small calculation leads to 



Solving for T, and comparing with the numerical values presented in Remark 12. 1[ 
gives 

Lemma 3.1. The minimal and maximal value for T such that the Harborth con- 
figuration exists are T — 0, and 



respectively. 

Using basic trigonometry, from this we are able to further deduce exact values 
for the crucial angle ip for extremal T: 

In case of T = 0, the points A, B, and C lie on one line. So do the points D, E, 
and F. Moreover the intersection point Z of these two lines together with C and D 
form an equilateral triangle, the sides of which have length 3. Angle f3 (see Figure 
2) becomes one of the angles of this triangle, and thus is equal to 60°. Furthermore 
the points D, F, and G form an isosceles triangle, with the length of the base side 
FG being one, and the other length being two. An analysis of the triangle formed 
by the lines DH, and the prolongations of DF and HJ, which contains the triangle 
DFG completely, allows us to calculate a. Since ip = a+(3 some further calculations 
show 

Corollary 3.2. ForT = 0, the angle ip in the Harborth configuration is the unique 
solution of 



in the interval [0,90°], which up to an error o/ lO^i'^ is 85.884964999269942°. 



(3.13) 



64r'^ - 567^ + 1 = 0. 



(3.14) 




(3.15) 
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Figure 3. Upper part of the Harborth configuration. 



As we have already observed, when T attains its maximal value, the points A, 
E, and F lie on one line, and form the side of the isosceles triangle ABE. Leaving 
the details to the reader, again only using basic trigonometry - and Mathematica 
for the calculation of trigonometric expressions - we are able to show 

Corollary 3.3. For T = ^\/7 — S-y/S, the angles a and (3 in the Harborth config- 
uration (see Figure 2) are the unique solutions of 

cos(/3) = ^ (^Vs + Vs - x/t-Vs^ , 

and 

68 + 3^230 + 34V5 + 9\/5 (s + ^230 + MV^) 
cos(a) — — — 

2(23 + 3\/5) y 97 - 3V5 + 3^230 + 34 V5 

in the interval [0,90°]. 

Since ip = a + f3, this leads to ip ~ 94.590425288952345° up to an error 0/10"^^. 

3.3. Prom the points D and F to the points H and J. Now we continue 
with our task of determining the minimal polynomial for that particular T, for 
which the Harborth graph exists, i.e., for which ip = 90° holds. Our trick is, not 
to calculate the coordinates of the points H and J directly in dependence of the 
second coordinate T of the point J5, but to introduce further variables X and Y 
which will lead to simpler equations. These new variables themselves will depend 
on the points D and F. 

For this, let us consider F as the origin of a new coordinate system, and the line 
ED as the new x-axis (see Figure 3). Let (s, S) be the coordinates of the point E 
with respect to this new coordinate system. Clearly, since DEE forms an isosceles 
triangle, D is described by the coordinates (2s, 0). Proceeding as above by using 
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Equations (I2.12p - (I2.17|) in this new context, we successively gelQ: 

(3.16) = 3- As^ + 4:SXG, 

(3.17) = 9 - 40s2 + iGs* + 16s2y2,, 

(3.18) = 9-48s^ + A8s^+{l2s-A8s^)xH + lGs'^x'jj, 
= -81 - 144s2 _ 352s4 _ 256s^ - 256s* + 

(3.19) (I44s2 + 896s'* + 2563^) yj^ ~ 256s^y%, 

(3.20) = 9- 36s2 + 16s'^+(l2s-16s3)xj + 16s^x}, 
= 81 - 504s2 + 1072s'* - 896s*' + 256s* + 

(3.21) (-144s2 + 256s'* - 256s*') yj + 2b6s%j. 

All these equations can be easily solved for the respective coordinates. In each 
case only one of the solutions is in accordance with our geometric construction. 
Below we present explicit formulas only for the coordinates of the points H and J: 



(3.22) XH 

(3.23) yn 

(3.24) xj 

(3.25) yj 



-3 + 12s2 + V3V-9 + 40s2 - 16s4 
3V3 + 4^/3s2 + V-9 + 40s2 - 16s4 
-3 + 4s2 - \/3\/-9 + 40s2 - 16s4 
-3V3 + 4%/3s2 + x/-9 + 40s2 - les^ 



»s 

Therefore the slope of the line HJ with regard to DF as a;-axis is 

(3.26) ma ■— = _ , . 

xj-XH 4s2 + V3V-9 + 40s2 - 16s4 

Now, again we consider the whole Harborth configuration: let new variables X 
and Y be defined by X := xd — xp and Y :^ yjj — yp, where now xd, Vd, and 
xp, yp are interpreted as the coordinates of the points D and F with regard to the 
initial coordinate system. Then the squared length of the line segment DF is given 
by X^ + Y'^. It follows that 

(3.27) 4s2=X2+y^ 
and Equation (|3.26p becomes 

3\/3 



(3.28) 



X2 



+ r2 + + 10 (X2 + _ (X2 + y2)2 



In order to be able to complete the Harborth configuration to the whole Harborth 
graph, the angle (f between the lines HJ and AC must be a right angle. On the other 
hand, we have = a + /3, where a is the angle between HJ and DF^ and (3 denotes 
the angle between DF and AC, as shown in Figure 2. Thus the equality a — 90° — /? 
must hold. Since 0° < a, (3 < 90°, we have tana = tan(90° - /3) = t^^^^- Thus the 



We advise the reader to keep in mind that, although we use the same notation as in Section 
12.21 now the coordinates have to be interpreted within the new coordinate frame. 
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respective slopes satisfy — l/mp. The slope mp of f3 is given within the first 
coordinate system by mp = Y/X. This, together with (|3.28|) . implies 

(3.29) 3V3 ^ X 

X^ + Y^ + ^-9 + 10 (X2 + r2) _ + y2)2 Y' 

which some calculations show to be equivalent to 

(3.30) 0= {X^ + Y^) (-27 + 30^2 - + gVsXY - iX^Y^ 



This implies 

(3.31) = 27 - 30X^ + iX* ~ eVsXY + iX^Y'^ . 
We set F{X, Y) 27 - 30^^ + iX'^ - 6\/3Xr + 4X2y2_ 

Next we produce polynomials Px.t and Py,t which describe the connection between 
T and the new parameters X, and Y, respectively. To do that, this time we use 
Mathematica's resultant and factorization facilities, as described above. Starting 
with the polynomial X — xd + xf, and the polynomials Pxo-T and Pxp,T given by 
(|3.3p and (|3.12p . we are thus able to successively eliminate the variables xd and 
Xf, and get 

Px,T '■= 

108T2-684r''+1344r''-1344r«+V3 (-36T + 300T^ - 616r^ + 1024T^) X 

(3.32) ^ (g _ 2i3r2 + 496r* - 1216T'') + {36T - OGT^ + 256T'^) 
+ (-9 + 24r2 - 64T*) X* e Z[V3][X,T]. 

In an analogous manner we deduce a polynomial Py,t in '^\y,T\, which for the 
Harborth configuration describes the connection between these two variables: 

Py,t := 

81-405r2+900r''-1008T'^+448r* + (54r- 54r^ +456r^ ~ SIST'^) Y + 

(3.33) ^ 207r2 - 624r'' + hl&T'^) + (-18T + 48r^ - 128T-"^) V'^ + 
{9 - + G'iT^) Y^ . 

Finally we repeat this procedure with the polynomials F{X,Y), Px.t, and Py,t 
to eliminate variables X and T. I.e., first we let Mathematica calculate the resultant 
of F and Px.t with regard to X. We will not present the result here; let it be 
enough to state that the result is a polynomial of degree 8 in Y, degree 32 in T 
and total degree 32, which, but for a constant factor, cannot be factored further 
by Mathematica, even when considered over Z[a/3]. With the help of Mathematica, 
we are able to calculate the resultant of this polynomial and Py.T with respect to 
Y. This leaves us with a polynomial in the single variable T of order 156. Strangely 
enough, this final polynomial is reducible over Z[\/3], its factors (up to an integer 
constant) being (2T + ^/3), {2T - VS), (64T4 - 24T2 + 9)^ and three other integer 
polynomials of degree 22, 28 and 80, respectively. Here, we need only list the 
polynomial of degree 22, since this is the one which has the y-coordinate T of the 
point B of the Harborth graph as one of its real rootsQ: 



^In fact it is the positive root of smallest modulus. 
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Theorem 3.4. The minimal polynomial for the y-coordinate T of the vertex B of 
the Harborth graph is 

Pt -492075 + 52356780r2 - 1441635408r'' + 12222052416r^ - 60567699456r^ + 
189747007488r^'^ -417660420096T1V607025037312T" -655053815808r" + 
446118756352r^® - 422064422912T^'' + 437348466688r2^ 

Since this polynomial is primitive, and irreducible over Z, it is the minimal 
polynomial of ys = T. Thus we have achieved our desired first main result. The 
rest of this paper is concerned with the determination of the minimal polynomials 
of the other coordinates, and some of their properties. 

4. Minimal polynomials for the ^-coordinates 

4.1. The points D,E, and F. The method to calculate the minimal polynomial 
for each of the y-coordinates of the points D, E, F is very similiar to what we 
have done in the last subsection. Let P denote one of the points D, E, or F. 
Using Formulas p.ip . (|3.7p . and p. lip , respectively, which describe the connection 
between the y-coordinate yp and the parameter T by way of irreducible polynomials 
Pyp.T G '^[yp.T] of absolute degree 2,2 and 4, we let Mathematica calculate the 
resultant of Pyp^r and Pt from above with respect to T. The result is a polynomial 
in Z[j/p]. It is of degree 44 in case of yn and yE and of degree 88 for yp. The first 
two of these polynomials each split over Z into two irreducible factors of degree 22. 
The polynomial for yp splits into one component of degree 44 and - again - into 
two polynomials of degree 22. In each case one of these factors corresponds to the 
actual y-coordinate. 

Theorem 4.1. The minimal polynomials for yD,yE o,nd yp are 

Py^ := -2470693585135788 + 1679453964496051893y|, - 2462573171102886288^1, + 

18471479139293280482/1, - 888334179987132288j/|, + 302241307009227264yl,° - 
74768143621533696yo + 13516084620361728^" - 1721332250836992j/i,^ + 
139442448236544yl,** - 6126808596480?yB + 109337116672?y|,^ 

Pyj, := -387038865725307 + 255845547796716j/| - 1080696123714384y| + 

985178573370432y| + 290816529555456y| + 12294226404679682/^? - 
399291497201664i/i? - 2269538689351682/^ - 1459143164559362/f + 
84049703993344j/|f - 94620310568962/1° + 4373484666882/f , 

Pyp := -6156736033068 + 4132620043369020y| - 280695352024663472/1- + 

54174190167055116y| - 443212523555443202/1 + 16893977313239424j/],'' - 
3430375146685440y]^ + 7819648176291842/]^ - 165954075623424t/].'^ -I- 
164009307013122/],* - 579898179584y|!' + 273342791681/1?. 

4.2. The points G,H, and J. When determining the minimal polynomials of the 
y-coordinates of the points G, H, and J, we have to keep in mind that in Section 
13.31 we used a different coordinate system for determining their coordinates than 
for those of the points A to F. Thus first we will have to transform the former into 
y-coordinates within our original system, which was centered in A. Since in this 
section it is paramount not to confuse these systems, in the sequel we will denote 
coordinates with regard to the coordinate system centered in F with capital letters, 
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and those with regard to the original one centered in A with smaU letters. Thus, 
e.g., Equation (|3.19p now becomes 

= -81 - 144s2 - 352s* - 256s^ - 256s^ + (1445^ + 89Qs'^ + 256s^) - 256s*Y^. 

We will have to treat Equations (|3.16p - (|3.21[) accordingly. With this notation, 
our goal has become to calculate the characteristic polynomials for j/g, i/h and yj. 
Elementary analytic geometry tells us that the connection between "old" coordi- 
nates {xp,yp) and "new" coordinates {Xp,Yp) can be described by 



(4.1) 



for P G {G,H, J}, where X — xd — xp and Y = yn — yp, as defined in Section 
Thus, 



(4.2) 2s{yp - yp) = Xp -Y + Yp ■ X. 

Setting zp := yp — yp for the moment, successively using Equations (|3.16p - (|3.2ip . 
which described the connection between new coordinates of the points G, i?, J, and 
the parameter s, and the equality 4s^ ~ X^ + Y^, after some calculations in the 
usual manner, we are able to deduce irreducible polynomials in Z[V3][zp, X, Y] for 
P G {G,H, J}. Using Equations (|3.32|1 and (|3.33p we are further able to ehminate 
X and Y, and deduce polynomials in Z[zp,T] of total degree 188 for each of the 
points. Each of these splits again, leaving us with irreducible polynomials which 
are of degree 20 in T and total degree 20 for the points G, J and degree 24 in 
T and total degree 24 for the point H. Resubstituting yp — yp ioi zp and using 
p. lip to eliminate yp^ we are left with integer polynomials in variables yp and T 
of degree 112 for the points G, J and one of degree 128 for H, which this time split 
off irreducible polynomials of total degree 20 (for G and H) and 24 (for J) , respec- 
tively. Calculating the resultant of these polynomials and Pp^ thereby eliminating 
T, once more in each case we get polynomials of degree 176 for yc, yn, and yj, 
respectively. Each contains among others an irreducible factor of degree 22 - the 
minimal polynomial, which we were looking for. Therefore we have: 

Theorem 4.2. The minimal polynomials for the y- coordinates of the vertices G, H, 
and J of the Harborth graph are 



Pya -912811377667500 -I- 16117998953248125y| - 36709013218422600j/g + 
37940201286814800yS - 234638874818542082/S + 10021184125203456i/g 
3290335763447808yG + 888521341648896yG - 192809455583232^^'' + 
29839017902080j/cf - 2742026240000yG? + 1093371 16672yg*, 



Py„ := -12148787578527675 - 123412000423046805^^ - 441020584930952232y|f + 
273168911377174014y|f - 27343071784237320i/|f - 36671 16898760364j/iF + 
8230449869876162/1,^ - 32095868573376?/]^* - 47799851427842/lf + 
615643279360?/3| - 27098808320j/|F -h 427098112?/!,^ 

and 

Pyj -9964518750000 + 57027771 1828125y J - 17805529663875002/5 + 

849106838377800y5 + 644904447905880y5 - 1020482802548282/}° - 
56106534718368y,}2 90274337585282/" -h605520976896y:}'' - 103349145600j/:^* - 
2815229952yf 4- 42 70981122/f . 
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5. Minimal polynomials for the ^-coordinates 

As we initially announced, we want to give minimal polynomials for all the 
coordinates of the most important vertices of the Harborth graph, where the center 
of origin is supposed to be the center of symmetry of the whole graph, i.e., the 
point K in Figure 2, and the axes are the axes of symmetry of the Harborth graph. 
This means that while the y-coordinates remain the same, when we shift our origin 
from A to K, the x-coordinates (with respect to the coordinate system centered in 
A) have to be shifted by —xj, i.e., we have to set 

(5.1) a;?,™ ^xp^xj 

for all points P, where the coordinates on the right hand side denote coordinates 
with respect to the origin A. 

As a further difficulty, again we have to pay attention that up above the x- 
coordinates for the points G, H and J were given with respect to yet another, third 
coordinate system, which had F as its origin and was rotated when considered 
within the other two coordinate frames. 

In the sequel, to avoid misunderstandings, we will switch notation, and denote all 
^-coordinates with respect to the system centered in A by Xp"^, those with respect 
to the system centered in K will now become xp. 

5.1. Coordinate transformations for the points A to F. Proceeding step by 
step as in Section 14.21 but starting from the second equation resulting from (j4.ip . 
i.e., 

(5.2) 2s{x°j'^ - xf'^) = Xj ■ X -Yj - Y, 

we are able to deduce an irreducible polynomial P^oU in Z[V3][x°}'^,T] of total 

degree 24, and finally the characteristic polynomial of x°j'^ of degree 22. Since it is 
an even polynomial, it describes the coordinate xa = —x'j'^ as well. Thus, 

Theorem 5.1. The minimal polynomial for the x-coordinate of the vertex A of the 
Harborth graph, where the coordinate system is the one given in Figure 1, is 

P^^ ■= -830376562500 -H358127000000a;i -34144387143750a;i -F96857243056800a;i - 
68697978132015a;5i - 189712941147x\'' -f 6188723588664a;Jf -704220643376a;Ji' - 
52577813248a;Jf + 27196394496xjf - 2918612992a;i° + 106774528a;f . 

With the help of the polynomial P ou rp from above, we can produce a poly- 

nomial Ptrafo G 'L[xp,x°p'^ ,T] that describes the coordinate transformation (|5.ip . 
by calculating the resultant of P^oid ^ and the polynomial xp — x°p'^ + x°j'^ with 

respect to x°j'^. This is of total degree 24, of degree 24 in T, and of degree 8 in both 
variables xp and x°p'^ . 

Applying the same method as above, i.e., first calculating the resultant of this 
"transformation polynomial" and the respective polynomials Pxp,t, which we now 
have to interpret as polynomials in the variables Xp'^ and T, followed by a factor- 
ization over Z[-\/3], and finally repeating this process with the resulting polynomial 
and Pt, we are able to deduce the characteristic polynomials of the x-coordinates 
xb to Xp- Since this procedure should be standard to the reader by now, we will 



^We call a polynomial even if it only contains monomials of even degree. 
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not go into details anymore, but will only present the results. Once again we stress 
the fact that these polynomials are for x-coordinates with respect to K as origin: 

Theorem 5.2. The minimal polynomials for the x- coordinates of the vertices of 
the Harborth graph, where the coordinate system is the one given in Figure 1, are 

P^^ -17372788157292129 +85946816541669534a;| - 172967171143553289x1 + 
125428630440736260x1 - 35361034276033728x1 + 4402034757921792xM' - 
436015591392256xf + 77220067192832x^* - 11054716223488x^® + 
874491412480xf - 34734080000xf + 557842432xf , 

P^^ -55268097000787592100 + 83653148035178006805xc - 
49933201015710366166x^ + 15170804748275250138xS - 

2623723693990622868x1; + 292733387369474292x^° - 24051159678783648x^^ + 
1563610131071808x^'' - 77064294460416x^^ + 2572257472512xjf - 
50083921920x^^ +427098112x^2^ 

P^^ := -15937557042969 + 69169635141939x13 - 133600085051911x1, + 

150590940104181x1, - 109441808559384x1, + 53597367271968xi? - 
17996039805696x1? + 4144963934208x1,* - 647005151232xif + 
66726690816xif - 4293132288x1? + 139460608x1?, 

P^^ ■= -30534686672400 + 184473995962680x| - 493600710483009x| + 

800738068318020x1 - 883225203916608x| + 687262746783744xM' - 
378024688788480x^2 ^ I45061641105408xj;' - 37695035736064xf + 
6218402758656x^* - 582162055168x|' + 27334279168xf , 

P^^ := -622521 + 20028276x|, - 150285424x|, - 349270464x|, + 7694997504x|, - 

5213620224x1? - 109200064512x1? + 709185896448x1? - 1112735219712x1? + 
387346071552x1? + 124822487040x|? + 8925478912x|?. 

5.2. The "Coup de grace" - the minimal polynomial for xo- Due to the 

complexity of expressions we were not able to use the above procedure to calculate 
the characteristic polynomial for xg - even with the help of Mathematica. Thus 
we have to resort to one final trick - yet another coordinate system. For this, first 
we mirror the Harborth configuration. After that we choose the point J as the 
new origin, and the ray JH as the positive part of the new a;-axis (see Figure 4). 
Consequently, the new y-coordinates are the x-coordinates from above. 

Let (m, U) be the coordinates of the point G with respect to this coordinate 
system. Using Equations ()2.1ip and (|2.12|) . adapted to the new coordinate frame, 
we get 

(5.3) - 3 + - AUyp + ^yl = 0, 

where yp now denotes the y-coordinate of F with regard to this new system. Since 
yp is equal to xp of old, we thus have calculated a polynomial Pxp,u- Calculating 
the resultant of this polynomial and P^p with respect to xp leads to a polynomial 
in Z[U] of degree 44. Factoring this, we achieve an irreducible polynomial of degree 
22, which is the minimal polynomial of U and thus of xq. 

Theorem 5.3. The minimal polynomial of the x-coordinate of the vertex G of the 
Harborth graph is 

P^c := -106929+9380331x^ -257190919xg +2410771629x^ - 11872837680x|. + 
35430882432x^1' -66974055936xJ;2 +79549160448xJ^ -56180293632xJ| + 
20514865152xJ| - 2573205504xg' + 139460608xg^ . 
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FiGURE 4. Part of the mirrored Harborth graph for the determi- 
nation of xg- Note that now H and J he on the x-axis 



Thus we have deduced minimal polynomials for all coordinates of the initial 
Harborth configuration, and consequently, because of the twofold symmetry and 
our particular, final choice of coordinates, for nearly all the vertices of the Harborth 
graph. This finishes our initial task. 

6. Coda 

With all the minimal polynomials at our disposal, we can finish this paper with 
some observations about their algebraic structure, and consequences for the Har- 
borth graph. 

First of all, starting with the minimal polynomials, using a computer algebra 
system we can once again show the existence of the Harborth graph - this time by 
algebraic means only: 

Theorem 6.1. For P g {A, . . . , J} let the coordinates {±xp, ±j/p) of points in the 
Euclidean plane be particularly chosen roots of the irreducible polynomials P^p , Pyp 
which were detailed in the previous sections. Then these coordinates satisfy the 
defining equations of coordinates of vertices of the Harborth graph. 

Proof. Since the actual calculations have to be done by the computer algebra system 
and cannot be presented here, the proof will just consist of a series of comments: 
First, with {ztxp,ztyp) we denote all four possible combinations of signs of the 
coordinates, taking into account the twofold symmetry of the Harborth graph, as 
well as the fact that the polynomials Pxp and Pyp are even. 

Second: Clearly Equations (|2.5p - (|2.18p have to be restated in accordance with 
our finally chosen coordinate system, which had the point K as its center. Thus, 
e.g., instead of (|2.5p one has to show that the coordinates of the vertices A and B 
in the first quadrant satisfy {xb — xa)^ + y\ = \^ i.e., the points are at distance 
1 from each other. Numeric approximations which determine our choice of roots 
have to be adapted as well. 
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Moreover calculations with the chosen coordinates are done in the sense of |9j, 
whereat the computer algebra system has to take the full brunt of the work. That 
is: since there are no explicit formulas for the particular roots of the polynomials 
PxpjPyp, which we use as coordinates, we have to understand these roots as com- 
pletely defined by their representing polynomials, and an isolating interval for each. 
Instead of an isolating interval one can equivalently use a numeric approximation 
of sufficient precision. In fact, this seems to be what the computer algebra system 
Mathematica does: it renders possible computations with algebraic numbers by way 
of root objects |15^ There, choosing a particular root actually means choosing one 
of the root objects produced by the Solve routine. This we did in accordance with 
the previously attained numerical results. The "proof" of the above theorem has 
been facilitated then by use of Mathematica's inbuilt RootReduce routine, which 
had to be applied to the defining equations of the Harborth configuration. In order 
to achieve that calculations are finished within an acceptable time, the equations 
have to be expanded and written as sums of products. More concretely, in the ex- 
ample above, after having assigned particular root objects to the variables xA, xB, 
and yB with the help of the respective polynomials, one successively has to calculate 
RootReduce [xA~2] , RootReduce [-2*xA*xB] , etc. until one is able to put everything 
together, thus to be able to calculate RootReduce [xA~2-2*xA*xB+xB~2+yB~2] , 
which indeed results in 1. □ 

Even though we have concentrated on one particular set of coordinates, our cal- 
culations have shown that, e.g., the y-coordinate ys of the vertex B in any case 
must satisfy a polynomial equation. Thus there are only finitely many possibili- 
ties for coordinates of B. Since by the geometric construction, which we used, all 
other vertices can be shown to depend uniquely on this one initial vertex and its 
embedding in the Euclidean plane, different embeddings of the Harborth graph in 
the plane - if they were to exist at all - cannot be transformed into each other in 
a continuous way. In other words: 

Theorem 6.2. The Harborth graph is rigid. 

Finally, closer scrutiny of the minimal polynomials leads to 

Lemma 6.3. Let zp he one of the coordinates of a vertex P of the Harborth graph 
different from zero. Then its minimal polynomial Pzp G is an even polynomial 

of degree 22 and signature (6, 8). That is, it has 6 real zeros and 8 distinct pairs of 
conjugate complex zeros. Conseguently we have, 

(6.1) P.p ^F.^oG 

where G = G '^[X], o denotes composition and F^p is a uniquely determined 
irreducible integer polynomial of degree 11 and signature (3,4). 

A precise proof for this lemma can be done by using any one of the existing 
algorithms for root isolation and some basic calculus. 

Theorem 6.4. The coordinates of vertices of the Harborth graph which are different 
from zero cannot be expressed in terms of radicals. Furthermore the Harborth graph 
as a whole cannot be constructed by compass and ruler alone. 
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Proof. A well known corollary in Galois theorjO tells us, that the zeros of a poly- 
nomial of odd prime degree, which is irreducible over a real number field, are ex- 
pressible in terms of radicals, if and only if either the polynomial has only one real 
zero, or all of its zeros are real. We have already ascertained that the polynomials 
Fzp appearing in Equation (|6.ip each have 3 distinct real roots. Thus the equations 
Fzp = 0, and consequently P^p = are not soluble by radicals. This proves the 
first assumption. 

We show the second assumption by an indirect proof: suppose that the Harborth 
graph were constructible by compass and ruler. This would imply that the vertex A, 
thus its x-coordinate xa and furthermore x\ would be constructible by compass and 
ruler. As a consequence of this, the order of the Galois group of the corresponding 
minimal polynomial, which is F^^, would be a power of 2, and the Galois group 
would be solublcQ. This, as we have just seen above, is absurd, completing the proof 
of the theorem. □ 
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See the corollary to Theorem 46 in [l]. 
'See Theorem 47 in [l]. 
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